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ABSTRACT 



We argue that the narrow Hne regions of Seyfert galaxies are powered by the 
transport of energy and momentum by the radio-emitting jets. This imphes that the 
ratio of the radio power to jet energy flux is much smaller than is usually assumed 
for radio galaxies. This can be partially attributed to the smaller ages (~ 10^ yrs) of 
Seyferts compared to radio galaxies but one also requires that either the magnetic 
energy density is more than an order of magnitude below the equipartition value, or 
more likely, that the internal energy densities of Seyfert jets are dominated by thermal 
plasma, as distinct from the situation in radio galaxy jets where the jet plasma is 
generally taken to be nonthermally dominated. If the internal energy densities of 
Seyfert jets are initially dominated by relativistic plasma, then an analysis of the 
data on jets in five Seyfert galaxies shows that all but one of these would have mildly 
relativistic jet velocities near 100 pc in order to power the respective narrow-line 
regions. However, observations of jet-cloud interactions in the NLR provide additional 
information on jet velocities and composition via the momentum budget. Our analysis 
of a jet-cloud interaction in NGC 1068, 24 pc from the core, implies a shocked jet 
pressure much larger than the minimum pressure of the radio knot, a velocity (probably 
accurate to within a factor of a few) ~ 0.06 c (18,000 kms~^) and a temperature 
of thermal gas in the jet ^ 10^ K implying mildly relativistic electrons but thermal 
protons. The estimated jet velocity is proportional to the jet energy flux and provides 
an independent argument that the energy flux in the northern NGC 1068 jet is much 
greater than previously supposed and is capable of providing significant energy input 
to the narrow line region. The jet mass flux at this point ^ 0.5Mq yr~^, is an order 
of magnitude higher than the mass accretion rate ~ 0.05 Mq yr~^ estimated from the 
bolometric luminosity of the nucleus, strongly indicating entrainment into the jet and 
accompanying deceleration. Consequently, the jet velocity near the black hole is likely 
to be mildly relativistic. We estimate an initial jet mass flux ~ 0.02 Mq yr^^ which is 
comparable to the mass accretion rate. This mass flux is consistent with the densities 
inferred for accretion disc coronae from high energy observations, together with an 
initially mildly relativistic velocity and an initial jet radius of order 10 gravitational 
radii. 
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1. Introduction 

The dominant paradigm for the excitation of emission hnes in the narrow hne regions (NLRs) 
of Seyfert galaxies invokes photoionization by a central power-law continuum. The apparent lack 
of ionizing photons in Seyfert 2s has been attributed to obscuration by a thick torus surrounding 



the central region of the Active Galactic Nucleus (AGN). Nevertheless, Binette, Fosbury, fc Parker 



(1993 ) have argued that there is a deficit of ionizing photons in all AGN leading to some scepticism 
of central power-law photoionization. On the other hand, it has been recognized for some time 
( Binette, Dopita, Sz Tuohy 1985D that excitation by autoionizing shocks in various classes of AGN 



can give a reasonable fit to the observed spectra and improve the agreement between theory 
and observation for the temperature-sensitive [OIII]A4363 line. Notable individual successes of 
this approach include a self-consistent model for the excitation of the emission line filaments in 
Centaurus A ( Sutherland, Bicknell, fc Dopita 1993 ) and the excitation of the disk in the central 



100 pc of M87 ( popita et al. 1997] ). In the latter case, detection by HST of high excitation lines in 
the ultraviolet together with the excellent agreement between the observed and predicted spectra 



make a strong case for shock excitation. [Bicknell, Dopita, Sz O'Dea (1997| ) have also shown that 



the narrow-line emission from Gigahertz Peak Spectrum (GPS) and Compact Steep Spectrum 
(CSS) radio sources can plausibly be related to the excitation of a dense ISM by the bow shock 
surrounding the radio source. That model has the additional feature of providing a natural 
explanation for the peak in the radio spectrum of these sources. Koekemoer fc Bicknell (1997| ) 



have also made a strong case for the shock excitation of the extended emission-line region in the 
radio galaxy PKS2356-61. For Seyfert galaxies, Viegas fc Contini (1997) have also suggested that 



emission lines from shocked gas in cloud-cloud collisions constitute a fraction of the total emission 
line output from the narrow line region (NLR). In view of the increasing interest in autoionizing 
shocks, a grid of models of shock excited spectra has been calculated (Dopita fc Sutherland 199"6a; 



Dopita Sz Sutherland 1996bD . These provide a good description of the observed spectra. 



The correlation between radio power and [OIII]A5007 luminosity in Seyfert galaxies 
discovered by [Whittle (1985|) over twelve years ago, suggests an intimate connection between 



the radio-emitting and emission-line components. Indeed, where observations are made at high 
enough resolution, detailed correspondence between radio and optical features is often found, 
along with evidence of dynamical structure and sources of local ionization indicative of shocks 
([Whittle et al. (1988[ ), [Goodrich (199^ , [Capetti et al. (1996[) ). In the specific case of NGC 1068 



Capetti, Axon, k Machetto (1997[ ) have shown, using IIST/WFPC2 [OIII] and IIa+[NII] images 
overlaid on a VLA radio image of the radio source, that there is a close correspondence between 
the radio and optical structure, strongly supporting the view that the NLR of NGC 1068 is the 
result of interaction with the radio plasma. Another key observation is that the radio axes of 
Seyferts correspond remarkably to the axes of their ionization cones, again suggesting that the 



radio and optical emission are causally related ( Wilson &: Tsvetanov 199^ ). The view advocated 



here is that the Whittle radio-optical luminosity correlation together with the unequivocal 
morphological correlations between radio and line-emitting gas is strong circumstantial evidence 
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that the mechanical energy and momentum transported by the radio plasma is responsible for the 
excitation of at least some, and probably most, of the narrow-line region (NLR). An alternative 
interpretation of a similar correlation for radio galaxies was given by Rawlings fc Saunders (1991| ) 



who suggested that the accretion disk luminosity (as indicated by the NLR [OIII] luminosity) and 
the mechanical luminosity of the jet are correlated. This suggestion was taken up by [Falcke &3i 
Biermann (1995| ) and Falcke, Malkan, fc Biermann (1995| ) who proposed that a substantial fraction 



of the energy dissipated by all AGN accretion disks finds its way into jet mechanical luminosity. 
However, this view does not really represent an alternative, since the inferred jet energy flux is 
comparable to the ionizing luminosity of the accretion disk and capable of having an important 
influence on the energy budget of the NLR gas. 

If the narrow line emission from Seyfert galaxies is to be powered by the expansion of the 
radio lobe, then, appealing to the GPS/CSS model of Bicknell, Dopita, fc O'Dea (199?] ) as a basis 



for order of magnitude estimates, the [OIII]A5007 luminosity is given by 

L{[OIII]) « LO X 10^2 f^-J\ Fe ergs s"^ (1-1) 



where 6 is the power-law index of the ambient density distribution and Fe is the combined energy 
flux in both jets. Hence, a median Seyfert [OIII] luminosity ~ ergs requires an energy 
flux in each jet ~ 5 x 10^^ ergs s~^. Thus, a major argument against the nonthermal plasma 
powering the emission- line region (e.g.. [Wilson (1997] )) is based upon the low radio luminosity of 
Seyferts (typically ~ 10^^ W Hz~^ at 5 GHz). Using standard conversion factors relating radio 
luminosity and mechanical jet power used for radio galaxies, one finds that the jet power falls 



short of what is required by at least two orders of magnitude. In earlier papers (e.g.. Bicknell 



(1995 ), Bicknell, Dopita, fc O'Dea (1997 )), we have represented this conversion factor by k^, the 
ratio of monochromatic power at frequency to jet energy flux. The advantage of using this 
parameter is that it is easily estimated using straightforward physics. The correlation between 
[OIII]A5007 luminosity and 1.4 GHz radio power is represented in Figure || where the correlation 



is plotted for a number of radio loud objects (used by Bicknell, Dopita, fc O'Dea (1997 )) together 



with the Seyferts in Whittle's sample. The lines represent the theoretical relationship between 



[OIII] luminosity and radio power predicted by the GPS/CSS model of Bicknell, Dopita, &: O'Dea 



1(1997 ). The values of logKi.4 (where kia represents Ki, at a frequency of 1.4 GHz) range from 



-14 (extreme left) to -10 (extreme right). Seyferts are separated quite clearly from the radio loud 
objects by, on average, three or even four orders of magnitude in K1.4 

In this paper, we examine the question of whether Seyferts are in fact characterized by such 
low values of ^1,4 and additional aspects of the energy and momentum budgets which have a 
strong bearing on the question of whether jets can power the NLR of Seyfert galaxies. Some of the 
following analysis is specific to NGC 1068 and the case of this galaxy, we examine the dyanmical 
link between the jet parameters and the expected parameters of the corona above the accretion 
disk in the vicinity of the black hole. 
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2. The Relation between Jet Power and Radio Luminosity. 

Traditionally, an effective value of Ky has been based on an assumed "efficiency of conversion" 
of jet power into total radio luminosity. However, this "efficiency factor" has a physical basis 
(iBicknell (19841) , |Eilek fc Shore (19891 ), [Bicknell (1995|) , [Begelman (199^) , [Bicknell, Do"^ite^ 



fc O'Dea (1997| )). For convenience, we summarize here the physics of this relationship as 



given by Bicknell, Dopita, fc O'Dea (1997 ). Consider a jet-fed radio lobe: The total energy, 
£^L accumulating in the lobe can be represented approximately as /ad -^e t where Fe is the 
jet energy flux, t is the age of the lobe and /ad ~ 0.5 is a factor which accounts for the 
energy lost adiabatically by expansion^. The synchrotron emission from the lobe is given by 
Pv ~ Csjn{a)K 3'^''+^'^/'^ zy-(«-i)/2^ y where Csyn is expressed in terms of the pacholczyck (197O0 



synchrotron parameters, q by Csyn(o) = 47rc5(a)c9(a)(2ci)("~-'^)/^, the electron energy density per 
unit energy, N(E), is given by N(E) = KE^"^, B is the (randomly oriented) magnetic field, v is 
the observing frequency, and V is the volume. The parameter K can be expressed in terms of the 
electron/positron energy density, eg via K = {a — 2)ee(7o?TieC^)^"~^^ [1 — (71/70)"^*^"^^]"^ where 
70 and 71 are respectively the lower and upper Lorentz factor cutoffs in the electron distribution. 
Writing the total electron/positron energy in the lobe as /e-E/L where /e is the electron/positron 
fraction of the total energy, we obtain for the ratio of monochromatic radio power to jet energy 
flux: 

K,«(a-2)Csyn(a)(7omec2)(«-2) \i - /e /ad i^^'^+'^Z^ t (2-1) 

The most important factors in this equation are (1) the electron/positron fraction, /e, (2) the 
magnetic field, B and (3) the age, t. This relationship is plotted in Figure |2| for a range of values 
of the evolutionary parameter r = fcfadt, assuming 70 = 1, 71 = 00, a spectral index of 0.7 
and an observing frequency of 1.4 GHz. The band of ^1,4 (—15 ^ logKi.4 ^ —13) which we 
suggest is relevant to Seyferts is indicated. For example a value of K1.4 10^^'^ and radio power 
~ 10^^ W Hz~^ indicates a combined jet energy fiux ~ 10^^ ergs s~^. From the plot it is evident 
that the values of B and r which are most conducive to K1.4 lying in this range are in the vicinity 
of -B ~ 10~^ — 10^^ G and corresponding values of r ~ 10^ — 10^ yrs. Typical equipartition 
magnetic fields and dynamical ages for Seyfert galaxies are approximately 10~^ G and 10^ yrs 
( Wilson 1997D , respectively, so that values of B and r in the required range are feasible. However, 



if r ~ 10^ yrs the magnetic field is required to be approximately an order of magnitude below 
equipartition. If the magnetic field is at equipartition we require r ~ 10^ yrs. The latter would 
be the case if r is less than the dynamical age as a result of dominance of the thermal pressure 
(compared to the nonthermal) making /e ^ 0.1. Thus this analysis raises the prospect that Seyfert 
jets, on the kpc scale, are dominated by thermal gas as distinct from radio galaxy and quasar jets 
which are normally taken to be dominated by nonthermal plasma. 

A related way of looking at this question is to ask whether the typical equipartition pressures 
inferred for the lobes of Seyfert galaxies are capable of driving the ~ 500kms~^ shocks necessary 



^/ad = 1 if there is no energy lost by expansion. 
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to explain the optical and UV emission line spectrum. The pressure required to drive a shock of 
velocity Vgh into gas with prc-shock density pps is PpsV^ ^ 6x 10^^ (nps/lOcm"^) (Vsh/500kms^^)^ 
where ripg is the prc-shock Hydrogen density. In order to produce the typical [SII] densities 
~ 10^~^ cm~^ inferred for NLR clouds we require rips ^ 10 cm~^ requiring driving pressures 
^ 10~^ dyn cm~^ well above the typical values of ~ 10~^ dyn cm~^ estimated for the lobes of 
typical Seyfert 2s from the nonthermal emission. Hence, we are again forced to the conclusion 
that either the lobes are either well out of equipartition or else that the lobe plasma contains a 
significant thermal component. 
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3. Implications for the ISM of Seyfert Galaxies 

Typically Seyfert radio sources are of the order of approximately a kpc in size. This statistic, 
combined with our inference of the energy flux in the jets required to power the narrow line region, 
has implications for the density of the confining ISM which we now calculate. These calculations 
are mainly concerned with the prototype Seyfert 2, NGC 1068, but they are easily applied to any 



other Seyfert. In this analysis we use the solution for the size of a radio lobe as given by Bicknell, 
Dopita, fc O'Dea (19971) (see also Begelman (1996 )). Take Xh to be the size of the lobe (from core 
to head), xq an arbitrary scale length, ambient density p = po{x / xo)~^ , -Fe the jet energy flux, ( 
the ratio of lobe pressure near the head to average lobe pressure and t the time, then 

Xh = xoe^^''''' (3-1) 

where 

(5-^ f^Etn_ (5-^)3 FE.43ti 

^ 187r(8-5) [poxlj • (8-6) no{xo/kpcf ^ ^ 

where 10"^^ -Fe,43 ergs s~^ is the jet energy fiux, Myr is the age, no cm~'^ is the ambient hydrogen 



density at a radius of xq, and, following Begelman (1996 ) we have taken C = 2. For a given size 
and age the implied density is 

^ (5-^)3C^ ^ (5 - Sf FE,43t^ (x,\ -(^-^) 

"° 187r{8-6) [poxlj ■ (8-,5) no(xo/kpc)5 Uo/ ^'^ 

For NGC 1068, the projected size is 700 pc and assuming that this is close to the actual size, we 
obtain Hydrogen densities of 3.1, 1.3 and 0.4 xFE,43ig cm~'^ for 5 = 0, 1 and 2 respectively. These 
densities are typical of those inferred for the disk ISM in the inner regions of spiral galaxies (e.g.. 
Wevers (1984 ), [Bosma (1981| )). That is, for any reasonable density profile, energy fiuxes in the 



range of 10'^^ — 10'^^ ergs s~^ and dynamical ages ~ 10^ yrs, we have an implied range of densities 
~ 1 — 10 cm^'^. This density is typical of what one would infer for most Seyferts and is in principle 
observable from the diffuse emission surrounding the lobes. Cecil, Bland, fc: Tully (1990| ) have 



determined densities 10^'^ cm~'^ from the [SII] lines emitted by ionized gas in NGC 1068. If 
these measurements represented an average density of the ISM in NGC 1068, then the implied 
jet energy fiux would be substantial (~ 10^^ ergss"^). However, the ionized gas is likely to be 
clumpy and if the gas is shocked as we suggest, then the [SII] density is likely to be about a factor 
of 50 — 100 higher than the pre-shock density, depending upon the magnetic parameter Bn^^^"^. 
This puts the pre-shock density into the range suggested by the above dynamical calculation. 
Thus, the existing ground-based measurements of the [SII] densities in NGC 1068 are consistent 
with our model but we cannot claim that they constitute definitive support for it. Nevertheless, 
this analysis and further calculations given below indicate that HST measurements of densities in 
this and other Seyfert galaxies are certainly of great interest. Mapping of the Faraday rotation 
structure may also lead to further insights. 
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4. Implications of the Energy Budget 



There have been numerous observations of radio emission in Seyferts which indicate colhmated 
outflow. However, only a few observations have the required sensitivity and resolution to actually 
detect jets. In the small number of cases that have been observed with sufficient resolution, e.g.. 
MKN 3 (iKukulaet al. 1993|) , MKN 6 ( [Kukula et al. 1996|) , MKN 78 (pedlar et al. (19891 ), [Whittle 
fc Wilson (1997| )), NGC 1068 ( [Gallimore et al. 1996| ; Pallimore, Baum, k O'Dea 1996|) and 
NGC 4151 (pedlar et al. 1993| ) the jets appear to be initially well-collimated and supersonic and 
then to disrupt and to become subsonic. In no cases do there appear lobes with bright hotspots 
as in classical FR2 radio galaxies. However, nearer to the core, jet knots occur, reminiscent of 
structures associated with shocks in supersonic jets. In NGC 1068, at least one of the radio knots 
near the core appears associated with an oblique shock where the jet is deflected by a narrow-line 
cloud. The northern NGC 1068 lobe may well be associated with the end of a subsonic jet 
and simulations relating to this situation would be of interest. Thus, dynamically Seyfert jets 
resemble FRl jets which appear to be initially supersonic and relativistic and to then undergo 
a transition to turbulent transonic flow. At this transition FRl jets are mildly relativistic with 
f3 = v/c ~ 0.6 — 0.7 ( Bicknell 1994 ; Bicknell 1995 ). The question then arises, are Seyfert jets also 
mildly relativistic in the region 10 — 100 pc from the core? 

Using the data in the above-cited papers it is possible to begin to address this question. In 
the following we assume that the [OIII] luminosity is the result of interaction of each jet with the 
ISM, i.e. the luminosity of the [OIII]A5007 line associated with one jet (i.e. one side of the NLR) 
is given by 

L([OIII]) « 0.02(1 -/ad) Fe (4-1) 

where /ad ~ 0.5 is the adiabatic factor mentioned above, the factor of 0.02 is the approximate 
fraction of the total emission line and continuum luminosity emitted in the A5007 line (cf. [Bicknell" 



Dopita, fc O'Dea (1997| )) and, as above, Fe is the jet energy flux . We use the energy flux inferred 



from this equation together with the [OIII]A5007 luminosities on the jet side to estimate the jet /? 



usmg 



3.6 X 10^2 



r-1 

p 

10~® dyn cm" 



10 pc 



r- 1 ■ 

1 + — F^X 



(4-2) 
(4-3) 



where p is the jet pressure, P = jet velocity/c, F is the jet Lorentz factor, ^jct is the jet 
cross-sectional area, and x = PcoidC^ /^P the ratio of cold matter energy density to relativistically 
dominated enthalpy ( Bicknell 1994] ). In order to obtain indicative estimates of jet velocity, we put 
p equal to the minimum pressure inferred from the radio images, taking the lower and upper cutoff 
Lorentz factors to be 10 and 10^ respectively^ We also take x = 1) since for an initially relativistic 



^The minimum pressure is relatively insensitive to these assumptions. The use of the minimum pressure here is 
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jet making the transition to subsonic flow, x ~ 1 ( Bicknell &: Begelman 1996| ). However, the 



estimated jet velocity is fairly insensitive to values of x near one. 

The indicative estimates of jet velocity and Lorentz factor implied by this procedure as well as 
the values of parameters used in the calculation are given in Table 1. The column labeled /([OIII]) 
in the table, is the fraction of the total [OIII] luminosity attributed to the jet. In all but galaxy 
the inferred jet velocities are of order 0.6 — 0.8 c. These calculations based on the energy budget 
alone, therefore could be taken to suggest that jet velocities in Seyferts are relativistic on the 10s 
of parsecs scale and that they are subsequently strongly affected by entrainment. Nevertheless, 
we discuss significant difficulties with such a proposition in the following section. These are 
principally based upon analysis of NGC 1068 radio and emission line data which suggest a much 
larger than minimum pressure in the jet and a larger momentum flux than could be provided by a 
relativistic jet. If this is the case for all of the jets that we have considered, then the real velocities 
are much lower. 

The exceptional case (Markarian 78) in Table 1 is of interest. Because of the distance of this 
galaxy, it is only possible to estimate the minimum pressure at ~ 270 pc from the core and our 
corresponding estimate of the velocity is 0.05 c, indicating, on the basis of the above assumptions, 
that the jet has decelerated to a subsonic velocity by this stage. The jet appears to be deflected by 
an [OIII]-emitting cloud at this distance from the nucleus ( [Whittle &: Wilson 1997 ) but there is no 



indication of a radio knot, indicating that the flow is subsonic, consistent with a smaller velocity. 
This deduction of deceleration is still valid if the pressure is dominated to the same extent by 
thermal plasma. 

As intimated above, energy considerations provide only a partial indication of jet velocities. 
Therefore, in the next section we consider the constraints placed on jet velocities by observations 
of jet-cloud interactions. This effectively brings the jet momentum budget into consideration. 



reasonable since, if the plasma is near equipartition, the total contribution of enthalpy and Poynting flux to the jet 



energy flux is similar to the contribution implied by equation 4-3 with p = pmin 
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5. The Momentum Budget from Jet-Cloud Interactions 

5.1. General considerations 

In a jet-cloud interaction as schematically indicated in Figure ^, a supersonic jet is deflected 
by a cloud producing an oblique shock in the jet. The post-shock pressure pressure drives a shock 
into the cloud which, if the cloud is dense enough, is fully radiative. The situation depicted in 
Figure f:jet-cloud is an oversimplification since the surface of the cloud will not remain flat but 
will be "gouged out" by the jet and also become filamentary - both as a result of the action 
of the Kelvin-Helmholtz instability. This will especially be the case if the post-shock flow is 
subsonic. In essence the force applied by the jet to the cloud is a fraction of the jet momentum 
flux and the estimates of shock luminosity presented below based on this idealized model probably 
underestimate this fraction. Nevertheless, the various parametric estimates we make below are 
probably indicative and motivate further numerical work on this topic. 

Given our idealized model, let us take p^h to be the post-oblique shock pressure in the jet 
and Pel to be the density of the cloud. The radiative shock driven into the cloud has a velocity, 
Kih ~ (4psh/3/Oc/)^''^ which together with the density and shock area determines the line fluxes 
from the cloud. Since Psh oc Pjct^et then the force exerted by the shock on the cloud, Psh^sh 
(where A^^^ cm^ is the shock area) is proportional to the jet momentum flux. We estimate the 
various coefficients of proportionality as follows. 

We parameterize the oblique shock by the pre-shock Mach number, Mjet, of the incident 
jet and the defiection angle, A^, of the jet, the latter being constrained by observations. We 
take 6i and 02 = 9i + AO to be the usual angles between the pre- and post-shock flow and the 
shock normal and Mx = Mjetcos^i to be the normal component of Mach number. Using the 
Rankine-Hugoniot conditions for a non-relativistic jet, one can readily show that 

2(M2-l)tan0i 

tanA6' = — ^— i^— (5-1) 

[(7 - 1 + (7 + 1) tan2 ^i] + 2 ^ ' 

where 7 = 5/3 is the ratio of specific heats. Given Mjet and A9 this equation can easily be solved 
numerically for 9i and then all of the other shock parameters, in particular, 



Psh i _^ 



cos^ 9i 



(5-2) 



The luminosity from shock excited emission lines is proportional to the shock area and we 
therefore need to relate the shock area to the area of the jet. As we show in Figure ^ this is 
determined by the addition of xi, the projection of the shock onto the cloud and X2, the distance 
traveled by the jet before a sound wave emitted from the edge of the jet intercepts the cloud. The 
total distance is given by 



xi +X2 = -Djet sec 9i (sin 6*2 -|- M2 cos 02 ) , 



(5-3) 
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where M2 is the post-shock Mach number, readily calculated form the Rankine-Hugoniot equations. 
The ratio of cloud shock area to jet cross-sectional area is given by j4sh/^jet ~ (xi + X2)/Djet- 
Expressing the jet momentum flux, approximately in terms of the energy flux by 



Pjet^et^jet 



2Fe 



jet 



1 -1 



1 + 



(7 - m?e 



jet. 



we have for the force exerted by the shock: 



where 



/(Mjet,A0) 



Psh 



jet 



1 + 



^et 

2 1 



jet 



sec 61 [sin 62 + M2 cos 62] 



and the jet velocity is given by: 



'jet 



jet 



2Fe 

Psh^shC 



/(Mjet,A0). 



(5-4) 



(5-5) 



(5-6) 



(5-7) 



The factor /(Mjet, A0) is easily estimated from the Rankine-Hugoniot conditions and a plot of 
this parameter as a function of Mach number for different deflections (suggested by the projected 
deflection at NLR-C) is given in Figure]^. Interestingly, the variation of /(Mjet, A^) with Mach 
number is minor and there is not a great variation with deflection angle. The physical parameter 
2Fe / (Psh^shc) can be determined from the energy flux required to power the optical emission 
excited by the lobe expansion and the parameters of the radiative shock in the dense cloud as 
shown below in the specific case of NGC 1068. 

In Figure ^ we also show the ratio of post-shock to pre-shock pressure as a function of Mach 
number. This is used in the following discussion of NGC 1068. 



5.2. Application to NGC 1068 

Using Merlin and VLA images together with an HST FOG [OIII]A5007 image of NGG 1068, 
Gallimore, Baum, fc O'Dea (199^ ) have shown that a significant number of emission line clouds lie 
adjacent to the Northern jet in NGC 1068, in particular clouds NLR-C, NLR-D and NLR-F. The 
cloud NLR-G lies along the projected direction of the jet which may be too faint at this position 
to show a knot. The projected areas of these clouds are all of order 10'*'^ cm^. These emission 
line clouds constitute an excellent example of where jet-cloud interactions are responsible for the 
excitation of at least part of the narrow line region and the following analysis aims to determines 
parameters for the radiative shocks which we assume exists in these clouds. In particular we are 
interested in the parameters of the cloud NLR-G near which the radio jet has a prominent radio 
knot (knot C). The jet also suffers a clear deflection ~ 23° at this location. Our interpretation of 
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the associated jet and narrow-line cloud morphology is that the deflection of the jet by NLR-C 
leads to the situation analyzed above, namely the production of an oblique jet shock and the 
driving of a radiative shock into the cloud by the associated pressure. 

For shocks in the velocity range of 500 — 1000 km s~^ the Ha luminosity is given in terms 
of the Hydrogen density nn cm~^, the shock area Ash cm^ and the shock velocity, V3 thousand 
km s~^, by 

L{Ha) = 5.3 X lO^^n^ A^t V^-*^ ergs s"^ (5-8) 

The pressure driving the shock, is given in terms of the cloud density, pc\ and the shock 
velocity, V^h; by p^h ^ Pci^h that the shock force is 

Psh^sh ~ PshVl « 4.5 X 10-6 Vf°-^^ L{Ha) (5-9) 

Here we use the observed velocity dispersion ~ 1000 km s^^ in the NLR of NGC 1068 as an 
estimate of the shock velocity. 

The Ha flux used in the above expression was estimated as follows: The HST Ha filter 
transmits (redshifted) [NII]A6548 and [NIIJA6583 and we took the transmission curve of this filter 
into account to obtain an estimated Ha flux of F(Ha -|- [A^//])/(l -1-0.74/2) where R is the ratio of 
[NIIJA6583 to Ha. The shock models in the velocity range 500 - 1000 kms-^ imply that R ^ 0.75 
so that we have used this value. Table 1 contains the measured Ha-|-[NII] and [OIII] fluxes in the 
brightest regions of the narrow line clouds as well as the [OIII]:Ha ratio. In principle one could 
use the velocity dependence of this ratio to estimate the shock velocity. However, the [OIII] flux 
is metallicity-dependent. Moreover, the high measured values of [OIII]:Ha suggest to us that the 
observed [OIII] line flux may be dominated by the shock precursor region. In table 3, therefore, we 
give estimates of number density, shock pressure and shock force for two shock velocities, 500 and 
1000 km s^^. The first two quantities are for a nominal shock area of lO^'' cm^. As equation ( ^-8| ) 
implies the number density estimate depends fairly strongly on the assumed shock velocity but is 
clearly in the vicinity of a few hundred cm~^ - typical of molecular cloud densities. The shock 
pressure is much less sensitive (see equation ( ^-9| ) and is in the vicinity of a few x 10"^ dynes cm~^ 
for shock areas ~ lO^'' cm~^. The shock force, of course, has the same velocity dependence as the 
pressure and is independent of the shock area. This is the parameter of primary importance here 
since it is the shock parameter which enters into the jet velocity estimate (see equation (|5-7| )). 

A clear prediction from these estimates is that densities measured from the [SII] lines should 
be in the vicinity of 10^^^ cm^'^ because of the density of the [SII] emitting region relative to the 
pre-shock density referred to earlier. 

Given the parameters in tables 1 and 2, the velocity of the jet near NLR-C is estimated to be 

2x10^"^ 

Act ~ g ^ ^Q34 ^ X /(Afjet, A0) « 0.1 X /(Mjet, A0) 

For a plausible range of deflections ~ 10° — 30° relevant to this particular jet-cloud interaction, 
/(Mjct, A6I) ~ 0.4 - 0.8 (see Figure ||) so that the jet velocity ~ 0.04 - 0.08 c. This velocity could 
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be significantly increased if we have underestimated the jet energy fiux required to power the 

NLR or if we have overestimated the force of the radiative shock responsible for the luminosity of 
NLR-C or if we have underestimated the fraction of jet momentum flux that drives the radiative 
shock. If, for example, the value of the parameter /(Mjet, A^) used above were a factor of two 
higher than we have estimated here, then the jet velocity would be a factor of two higher. Thus a 
jet velocity, vi 0.2 — 0.3 c, before knot C would not be out of the question. In view of this factor 
of two uncertainty in the estimation of physical parameters further numerical simulations of this 
process would be useful as would more detailed spectral observations of the NLR of NGC 1068. 
The latter would assist us to better estimate the parameters of the radiative shocks involved. 
Despite these reservations, this calculation demonstrates how the combination of total narrow- line 
region luminosities and individual cloud luminosities associated with jet-cloud interactions can be 
used to derive jet velocities. It will be interesting to see similar calculations applied to a sample 
of Seyferts in which the individual uncertainties in various parameters would statistically average 
out. We also add that it would be unlikely that we have estimated quantities so poorly that the 
above velocity is in error by an order of magnitude and that the jet velocity is in fact close to c 
and that the Lorentz factor is large. 

Another important feature arising from this analysis is that the shock driving pressure is 
larger, by about a factor of 40, than the minimum pressure inferred from the radio image and this 
factor is so large that even given the uncertainties of the above analysis, it is significant. This 
cannot be an effect of resolution since the shocked region of the jet adjacent to the cloud should 
be extended by 2-3 jet diameters in the jet direction (the observed elongation based upon the 
20 cm. image is about 2 jet diameters) so that the entire post-shock region should be resolved. 
Moreover, departures from minimum conditions of this amount are unlikely. Rather, this is 
another indication that there is an additional component of the jet pressure. We can now consider 
the implications of this by returning to the energy fiux. Constraints derived from this can be 
used to estimate the jet Mach number, as follows. For a nonrelativistic jet the energy flux can be 
expressed in the form: 



7 — i 



i + 2_1m.2 



2 



jet 



(5-10) 



7' 

where 7 is the ratio of specific heats. For a pressure, p = W~^p-Q dyn cm~^ Vjet ~ 0.06 c, 
Djet ~ 18 pc and 7 = 5/3, 

Fe « 1 X 10^3 + ^ j . 

The pressure here is the pre-shock pressure, for which we have no direct estimate. However, 
given that the post-shock pressure ~ 6 x 10~^ (^sh/lO'^'^cm^) ^ and that the jet energy flux, 
Fe ~ 10^^ ergss"^ the ratio of post-shock to pre-shock pressures, P2/P1 is given by 

where C ^ 0.6. The area of the radio knot adjacent to NLR-C (?a 5 x 10^^ cm^) indicates that the 



- 14 - 



actual area of the shock may be less than the nominal shock area cm ). Therefore, in the 
right panel of Figure ^ we have superposed plots of equation ( ^-12D ,for C = 0.6 and 1.2, on the 
plots of the theoretical shock pressure ratios. Allowing C to vary by a factor of 2 also allows for 
errors in the estimates other quantities such as the velocity and energy flux which are not known 
to better than a factor of 2. 



It is evident from the intersection of equation ( 5-12| ) with the calculated shock pressure ratios 
that a low, but supersonic jet Mach number ~ 1 — a few is favored. This is consistent with the 
morphology of this jet: The occurrence of knots, which it is natural to associate with shocks, is 
indicative of supersonic flow; however the jet spreads at the rate that we associate with turbulent 
transonic flow. Moreover, it is not as well-collimated as we would expect of a highly supersonic jet. 

Given the Mach number and velocity, the sound speed and hence the temperature of the jet 



(/zmpcV(7fc)Mj-2/32^ 



3, for example, T 



''Mrjpf^^ K. For 
1.6 X 10^ K. Thus the 



can be estimated, with the result that T 
/3jet ~ 0.06, T ^ 1.5 X 10^°Mj-2 K and if Mjet 
electrons are mildly relativistic and the protons are subrelativistic, assuming that they are at the 
same temperature. Thus, although the particle energies are not ultrarelativistic, they are extreme, 
and must be related to the plasma environment near the black hole. 



The mass flux, M, in the jet is also of interest. This is given by 

^1 



M 



2Fe 

y2 



1 + 



(7 



1)M2 



jet 



1 Mq yr- 



-1 



10^4 ergs s- 



1 + 



(5-13) 



For values of the jet Mach number, Mjet = 1, 2 and 3, the implied mass fluxes are 0.3, 0.6 and 
0.7 Mq yr~^, respectively. These estimates are all about an order of magnitude larger than the 
mass accretion rate rate, ~ 0.05 Mq yr~^, onto the black hole estimated from the bolometric 
luminosity ~ 3 x lO^'* ergs ( pallimore et al. (1996| ) corrected for our assumed distance of 
15 Mpc) and assuming an efficiency factor of 0.1. This accretion rate is close to the estimate 
0.04 Mq yr~^ of Maloney (private communication) derived from the properties of the maser 
emission. The discrepancy between the jet mass flux and the accretion rate is a clear indication 
that most of the jet mass flux is the result of entrainment and given the substantial amount of 
dense matter in the circumnuclear environment of NGC 1068, this is not surprising. This raises 
the question as to whether the initial velocity of the jet could be relativistic. To answer this we 
consider the implications of momentum and energy conservation in a relativistic entraining flow. 
In view of our above deductions concerning the composition of the jet plasma, we assume that the 
jet is thermally dominated and that initially, the energy and momentum fluxes are dominated by 
the rest-mass inertia. Hence the initial energy and momentum fluxes {F^.i and Fp^i respectively) 
are given by 



Fe,i ~ (ri-l)Mic2 



(5-14) 
(5-15) 



where /3i is the initial value of v/c, Ti is the initial bulk Lorentz factor and Mi is the initial flux 
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of rest mass. The energy flux is conserved, irrespective of entrainment. If the ambient pressure 
gradient does not have a major effect on the momentum flux (as is the case when the Mach 
number is high), then it also is approximately conserved. Thus the ratio of these two quantities 
is also conserved. When the flow becomes subrelativistic the ratio of energy to momentum flux is 
Vjet/2(l + 3/Mj2t) so that 

(Ti - l)c 1 / 3 



and the initial velocity is determined by the velocity and Mach number at knot C by 

Ti - 1 1 + 3/ML 

^^~/3jct (5-17) 



The nonrelativistic limit of this equation is 

These equations imply that for a transonic Mach number ~ 1 — 2 near knot C the initial velocity 
of the jet is about 2-4 times higher, i.e. about 0.1 — 0.3c, given our above estimates. As we have 
outlined above, it is not unreasonable to envisage "factor of order unity" changes in the parameters 
of the radiative shock, jet energy flux and the fraction of the jet momentum flux driving the 
radiative shock that have gone into the estimate of the jet velocity at knot C. Moreover, if the 
background pressure gradient increases the momentum flux, in the transonic regime of the jet, then 



equation (5-17) underestimates the jet velocity on the tens of parsecs scale, with the result that 



the initial jet velocity could be higher. (This can be seen more clearly from equation ( 5-16| ) with 



~ replaced by ;^.) Thus it is reasonable to conclude that the initial jet velocity in NGC 1068 is 
mildly relativistic and this is what we expect from a flow ejected from, say within 10 gravitational 
radii of a black hole. However, it is unlikely that the Lorentz factor is of order a few-10, as in 
radio galaxies and quasars. 

The following considerations of the mass flux enable us to relate the jet properties to that 
of the corona above the nuclear accretion disk. The initial mass flux, M ~ FEC~^(ri — 1)~^, 
i.e. 2F^E{c(3i)~'^ ?s 4 x 10~^/3~^Fe,44 Mq yr'^, in the non-relativistic limit. If the jet is mildly 
relativistic (say j3 ~ 0.5) then the implied mass flux ~ 0.02 Mq yr~^, close to, but less than, 
the above-estimated mass accretion rate ~ 0.05 yr^^. (The two are equal for j3i ~ 0.3.) We 
suppose that the jet originates from the coronal region of the black hole accretion disc and that 
the appropriate radius is of order 10 gravitational radii since this is where most of the accretion 
disc and related coronal dissipation occur. It is now reasonably well-established that the hard 
X-ray emission from Seyferts is thermally Comptonized emission from an accretion disc corona 



with an electron temperature, Te ~ 50 — 100 keV and Thomson optical depth, ~ 2 (Johnson 



et al. 199^ ). Adopting a coronal radius ~ 10 gravitational radii, then the electron density in the 
corona, 

nc ~ 1.0 X 10^^ rT(rc/10rg)-^Mf^ cm'^ (5-19) 
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where = GM/c^ is the gravitational radius for a black hole of mass M, and WMjMq is the 
black hole mass. Assuming that some of the corona is ejected in the form of a mildly relativistic 
jet of density rii, radius ri ~ lOr^, and a Lorentz factor Ti then the corresponding mass flux, 



Ml TTrfniTiPic 



2 



3.5 X 10^2 / ^ I / ^ I Mj (—] Tif3i Mq yr"^ (5-20) 



lOrg J \ Wvg J yric 

Given that Mj ~ 3 for NGC 1068, assuming that r ~ 2, that the jet velocity is mildly relativistic 
(again for argument, say /5i ~ 0.5) and also assuming that the jet density rii ~ ric then the 
initial jet mass flux ~ 0.1 Mq yr~^. This estimate is close to the mass flux estimated from the 
jet dynamics and any discrepancy could be attributed to the decrease of density from the corona 
to where the jet terminal velocity is established. Assuming that the jet and coronal radii are 
similar (i.e. ri ~ Vc), then our estimate is not very sensitive to our assumptions concerning 
these parameters. This good order of magnitude agreement motivates closer examination of the 
relationship between discs, their coronae and the jets ejected from them. 

The other issue to address here is the confinement of a jet with a pressure ^ dynes cm~^. 
This is too large a pressure to be provided by a hot interstellar medium; a number density 
^ 700 cm^'^ would be required and the cooling from this would be catastrophic. The other 
possibility, suggested by the observations of Cecil, Bland, fc Tully (1990| ), is that the jet is 



inertially confined by a dense cool outflowing wind. If the wind has a density ~ 100 cm ^, the 
overpressured jet drives shocks into it with a velocity 

, / P / ^ xl/2 

Vsh = 650 km s'^ —p-. =? (5-21) 



10 ^ dynes cm ^ / V 100 cm ^ 



There is evidence for such shocks from the observations of Axon et al. (1997 ) who find spectral 
evidence for line splitting of approximately 1000 km s~^ in the gas coincident with the jet (in 
projection). The presence of the coronal line of [FeVII]A3769 in their spectra is indicative of high 
excitation. Using the data of Pecil, Bland, fc Tully (199C| ) which indicate a velocity ~ 1000 km s~^ 
and an opening cone angle for the wind ~ 40°, the mass flux is Mwind ~ 0.8 Mq yr^^. (We 
have evaluated the cross-sectional area at 24 pc, the location of NLR-C.) Again, it is interesting 
that this mass flux is of order the mass accretion rate. The energy transported by such a wind 
~ 3 X 10^^ ergss"^ and therefore the wind is not competitive with the jet in exciting the narrow-line 
region. 

Finally, we also note here that the jet velocity implied by a conventional (radio galaxy) value 
of the parameter n^, would be unreasonably small. The value of ki 5 implied by our calculated jet 
energy flux ~ 10^^ ergs s~^ and the Wilson &: Ulvestad (199^ ) radio flux density, 3.8 Jy, for the 



linear part of the source is K1.5 ~ 5.5 x 10 Hz ^. If K1.5 were a factor of 10'^ — 10^ higher, say 
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10~^^ — 10~^^ Hz~^, the jet energy flux would be a factor of 10^ — 10^ lower and the jet velocity, 
implied by equation (5-7), would be in the range of 3 — 30 km s^^. Therefore, our interpretation 
of the jet and associated emission- line morphology, provides additional support for a much lower 
value of in NGC 1068 and consequently for much greater energy input into the NLR than 
previously supposed. The only way in which this conclusion could be in error would be if the cloud 
NLR-C were photoionized and its excitation had nothing to do with the deflection of the jet. 
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6. Discussion 

In this paper we have addressed the physics which is necessary in order to support the view 
that the NLR emission from Seyfert galaxies be powered, at least in part, by the energy and 
momentum associated with the radio jets. One requires jet energy fluxes ~ 10^^"^"^ ergs s~^ 
and a much smaller ratio of radio power to jet energy flux than is normally invoked for radio 
galaxies. A lower value of this ratio is suggested by the smaller ages of Seyferts compared to 
radio galaxies. However, we also require that either the magnetic fields in the lobes of Seyferts be 
sub-equipartition or more likely, that the fraction by energy of relativistic plasma in the lobes be 
much smaller than in radio galaxies. This raises the question as to where the composition of the 
jets is established. Are the jets initially relativistic in composition and subsequently diluted by 
entrainment of thermal plasma from the ISM or are they initially thermally dominated outflows 
with a small proportion of embedded relativistic gas? If we adopt the first point of view and 
assume that all AGN jets have initially relativistic velocities and that their internal energies are 
initially dominated by relativistic plasma, then we infer mildly relativistic velocities in four Seyfert 
jets mostly within 100 pc of the core. On the other hand, in NGC 1068, additional information 
provided by a jet-cloud interaction near the core, suggests a lower velocity ~ 0.04 — 0.0.08 c 24 pc 
from the core and, in order to maintain the same energy fiux, a substantial contribution to the 
jet pressure from thermal gas is required. As we noted above, the velocity estimate from the 
parameters of the jet-cloud interaction also support a significantly higher jet energy flux in the 
northern NGC 1068 jet than previously supposed and therefore support our suggestion of a much 
lower value of the parameter in Seyfert galaxies. A counter to this argument would involve the 
proposition that NLR-C is, in fact, photoionized, and that the deflection of the jet at this point 
has nothing to do with its excitation. This could be resolved by detailed HST optical and UV 
spectra of this cloud which could establish whether it is photoionized or shocked (see, for example, 
Dopita et al. (1997|) ). 

At present the notion that Seyfert jets may not be initially ultrarelativistic in composition, 
is only based upon the analysis of data on one galaxy, NGC 1068. However, the alternative, 
that Seyfert jets are dominated by thermal plasma from the outset, is probably more sustainable 
as a general proposition. If the internal energy of Seyfert jets were relativistically dominated, 
then one would have to explain why a ~ lO^'* ergs s~^ jet of relativistic plasma produced so little 
radio emission on kpc scales. Dilution by substantial entrainment of thermal gas and a related 
energy loss by the relativistic component as it mixes with the thermal plasma may go part way 
to explaining the low level of extended radio emission. However, a model of this sort may meet 
substantial physical difficulties. It is far more straightforward to invoke a substantial thermal 
component from the outset. Nevertheless, whilst we have shown, on the basis of our interpretation 
that the NGC 1068 jet is probably not dominated by ultrarelativistic particles, the inferred 
temperature of the NGC 1068 jet, at 24 pc from the core, ~ 10^ — 10^'^ K is substantial and must 
be related to the environment of the black hole where it originated. As we have shown, a high 
pressure jet must interact substantially with its surroundings via shocks and the high excitation 
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line emission coincident with the jet, discovered by Axon et al. (1997 ), is good evidence for this. 
The mass flux in the wind, ~ 1 Mq yr~^ is similar to that in the jet at that distance from the core, 
and is similar to the mass accretion rate onto the black hole. 



Work by Colbert (1997) supports the view that Seyfert radio plasma is, at least on large scales. 



dominated by thermal gas. He finds that, in Seyferts with large-scale outflows, the overwhelming 
contribution to the lobe pressure is provided by soft X-ray emitting gas, the radio emitting plasma 
contributing about 1% consistent with our requirements on the energy flux. Dominance of the 
lobe pressure by (T ~ 10'' K), thermal gas has another physical consequence which enhances the 
self-consistency of our model: The isobaric cooling time for shocked gas, fcooi ~ S.GpZgT^'^^ Myr 
(where 10~^p-9 dyn cm~^ is the pressure and 10^ T7 K is the temperature) indicates lobe cooling 
time scales ~ a few Myr. This provides an explanation for why the dynamical ages of Seyferts 
appear to be ~ a few Myr (corresponding to characteristic expansion velocities ~ 500 km s~^ and 
sizes ~ kpc). Once the gas in the lobe cools, the lobe collapses and then builds up again over 
this timescale. Obviously, the gas in the initial region of the jet needs to cool from ~ 10^ K both 
adiabatically and through mixing in order to reach ~ 10^ K. 



The paper by Nelson &: Whittle (1995| ) is also of interest in this context. They have shown 



that when the total powers of Seyfert and radio galaxies are plotted against bulge magnitude, the 
radio galaxies lie well above the sequence defined by the Seyferts. However, restriction to the core 
powers of radio galaxies sees the radio galaxies continuing the Seyfert sequence. That is, when the 
radio powers of Seyfert and radio galaxies are compared on similar scales, they form part of the 
same sequence. This suggests that entrainment on the kpc scale may have some influence on the 
level of kpc-scale radio emission. However, the Seyferts do have, on average, core powers which 
are lower than those of the radio galaxies consistent with them being weaker radio sources from 
the outset. 

Since the mass flux estimated for the jet at 24 pc from the core exceeds the accretion rate 
into the black hole by about an order of magnitude, it is likely that the jet mass flux at this point 
is the result of entrainment. Our estimate the jet velocity near the black hole is subject to the 
uncertainties in the estimates of jet velocity and Mach number near knot C. However, it is not 
hard to justify a mildly relativistic velocity. The jet velocity near the black hole is also related to 
the initial mass flux. If the initial velocity is mildly relativistic then the mass flux is similar to but 
less than the accretion rate. Moreover, the mass flux is consistent with a typical Seyfert coronal 
density ~ 10^^ cm~^ being ejected at a mildly relativistic velocity. Thus consideration of the NLR 
excitation has led us to a fundamental linkage between jet properties and the coronal properties 
of accretion discs. 

The remaining point we wish to make is that, for Seyfert jets dominated by thermal plasma, 
one generally expects a large amount of internal Faraday rotation. For example at knot C in 
NGC 1068, we calculate a rotation measure of 20,000 radm^^ for an equipartition field ~ 10"'^ G. 
Unfortunately internal jet depolarization is diflicult to disentangle from depolarization by the 
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NLR. 
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Table 1: Jet velocities inferred from [OIII] and radio data 



Galaxy & 


Distance to 


Diameter 


Pmin 


fiioni]) 


L{[OIII]) 


/? 


Radio reference 


Component 


nucleus (pc) 


(pc) 


dyn cm~^ 




ergs s"^ 




Optical reference 


Mkn 3 E 


69 


41 


2.6 X 10"^ 


0.5 


4.8 X 10^^ 


0.63 


Kukula et al. (1993) 


Mkn 3 W 


141 


36 


2.0 X 10-^ 


0.5 


4.8 X 10^^ 


0.73 


Koski (1978) 


Mkn 6 4 


127 


51 


5.8 X 10-^ 


0.5 


1.7 X 10^2 


0.63 


Kukula et al. (1996) 
Koski (1978) 


Mkn 78 W 


271 


180 


5.7 X 10-^ 


1.0^ 


8.7 X 10^^ 


0.05 


Whittle k Wilson (1997) 
Whittle & Wilson (1997) 


NGC 1068 C 


24 


18 


1.4 X lO"'^ 


1.0^ 


1.0 X 10^2 


0.77 


Galhrnore et al. (1996) 
Storchi-Bergmann et al. (1995) 


NGC 4151 CI 


114 


37 


9.5 X 10-9 


0.5 


2.2 X 10^^ 


0.71 


Pedlar et al. (1993) 


NGC 4151 C2 


56 


20 


1.4 X 10-^ 


0.5 


2.2 X 10^^ 


0.73 


Anderson (1989) 


NGC 4151 C5 


114 


23 


1.4 X 10-^ 


0.5 


2.2 X 10^^ 


0.82 





^ [OIII] luminosity of western region of source. 

^ The entire [OIII] luminosity of NGC 1068 is assigned to the northern region because of the one-sided ionization cone. 



Table 2: Observed Ha and [OIII]A5007 luminosities 



Cloud F{Ha + [NII]) F{Ha) log L{Ha) F{[OIII]) log L{[OIII]) [OIII]/Ha 

ergs cm~^ ergs s~^ ergs cm~^ s~^ ergs s~^ 

NLR-C 8.50 X 10-^3 5.47 x IQ-^^ 40.14 2.03 x lO'^^ 40.71 37 

NLR-D 4.01 X 10-^3 2.58 x lO^^^ 39 gi 9.00 x 10"^^ 40.35 3.5 

NLR-F 4.40 X 10"^^ 2.83 x lO"^'^ 39.85 1.50 x lO'^^ 4Q 57 5 2 

NLR-G 2.59 x 10"^^ 1.67 x lO'^^ 39^52 9.99 x 10"^^ 40.40 6.0 
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Table 3: Estimated number densities and shock pressures 

Cloud V^i, nH(A,h/1040 cm2) psh^sh Pshl^sh/lO^o cm^) 

kms~^ cm~^ dynes dynes cm~^ 

NLR-C 500 1500 8.2 x 10^^ 8.2 x 10-^ 

NLR-D 500 690 3.8 x lO^'' 3.8 x 10"^ 

NLR-F 500 750 4.2 x 10^^ 4.2 x 10"^ 

NLR-G 500 440 2.5 x 10^^ 2.5 x 10"^ 



NLR-C 1000 280 

NLR-D 1000 130 

NLR-F 1000 140 

NLR-G 1000 80 



6.2 x 10^4 6.2 X 10-6 

2.9 X lO^'' 2.9 X 10-6 

3.2 X 10^^ 3.2 X 10-6 

1.9 X lO^"^ 1.9 X 10-6 
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Figure captions 



Figure |^: Predicted [OIII\ luminosity as a function of 1.4 GHz radio power for values of 

logKi.4 = —10, —11, —12, —13 and -14 overlaid on data for radio-loud objects and Seyfert 
galaxies. Data are from Gelderman &: Whittle (1996D (GW), [Tadhunter et al. (1993| ) and 



Morganti, Killeen, fc Tadhunter (1993| ) (collectively referred to as TM) and [Whittle (19"85| ) 



Legend: Filled circles - CSS sources from GW; filled squares - CSS sources from TM; open 
circles - GW FR2 radio galaxies; open squares - TM FR2 radio galaxies; diagonal crosses - 
GW QSOs; plus signs: TM compact flat spectrum sources; open triangles: TM FRl radio 
galaxies; filled hexagons: Seyfert galaxies from Whittle (1985| ). Upper limits are indicated in 
the usual way. 

Figure |: The ratio Ki.4 of radio power at 1.4 GHz to jet energy flux for different values of the age 
parameter r = fefadt- A spectral index of 0.7 and a lower cutoff, 70 = 1.0 have been 
assumed. However, the value of Ku is not very sensitive to these parameters. 

Figure |: This fi gure defines the parameters of a jet-cloud interaction. The extent of the high pressure 
post-jet-shock region is determined by the sound wave which propagates towards the cloud 
from the edge of the jet. This high pressure region drives a radiative shock into the cloud. 

Figure ^: The left panel shows the factor /(Mjet,^i) defined by equation ( ^-6| ) which parameterizes 
the force of the post jet-shock region (psh^sh) in units of FE/Vjet for jet deflections (10°, 
20° and 30°) as indicated. The right panel shows the pressure ratio (solid lines) P2/P1 as 
a function of pre-shock jet Mach number for the given deflections. Superimposed on these 
curves are dashed curves {P2/P1 = C(l + Mjg^/3), C = 0.6 and 1.2) representing the pressure 
ratio implied by the radiative shock analysis and the jet energy budget for NGC 1068. Both 
panels are for a non-relativistic jet with the ratio of specific heats 7 = 5/3. 



-28- 




log Plot (W/Hz) 



Fig. 1.— 
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Fig. 2.— 




Fig. 3.— 
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Fig. 4.— 



